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Abstract. We present a thorough analysis of the electronic detection of charged 
particles, confined in a Penning trap, via image charges induced in the trap electrodes. 
Trapping of charged particles in an electrode structure leads to frequency shifts, which 
are due to image charge and space charge effects. These effects are of importance for 
Penning trap experiments which involve high charge densities or require high precision 
in the motional frequencies. Our analysis of image charges shows that only (higher 
order) odd powers of the particle displacement lead to induced charge differences, giving 
rise to a signal. This implies that, besides the centre-of-mass frequency of a trapped 
particle cloud, also higher order individual particle frequencies induce a signal, which 
can be picked up by an electronic detection circuit attached to the trap electrodes. 
We also derive analytic expressions for the image charge and space charge induced 
frequency shifts and perform simulations of space charge effects. In relation to this, 
we discuss the consequences of the shifted particle frequencies for resistive cooling of 
the particle motion. 
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1. Introduction 

Penning traps have proven to be versatile tools for extended investigations of localised 
particles under well defined conditions, thus allowing for high-precision experiments. 
Such experiments often rely on precise particle manipulation and require knowledge 
of the motional frequencies and other dynamic properties of the trapped particles. 
However, the presence of charged particles in a Penning trap is accompanied by image 
charge and space charge effects, which change the effective trapping potential and 
therefore lead to oscillation frequency changes. Image charge effects are caused by 
particle-induced image charges in the trap electrodes and in most cases lead to only 
small frequency shifts. However, in precision experiments they still need to be considered 
[TJ Space charge effects are caused by the charge density of trapped particles and 
lead to considerable frequency shifts and broadening [H El El IH El E] • 

In the following we discuss ideal particle motions in a Penning trap. We present a 
thorough analysis of the formation of image charges and determine image charge and 
space charge induced frequency shifts of the particle motions. Electronic detection of 
the axial particle frequency (section 5) will be discussed, followed by a section on the 
implications for resistive cooling of the (shifted) axial frequency (section 6). The aim 
of this study is to give a consistent quantitative treatment of these effects. Similar 
effects also occur in other types of traps, such as Paul or radio-frequency traps, but 
their analysis is beyond the scope of this study. 

2. Single particle motion in a Penning trap 

In a Penning trap, axial (z) confinement of a single particle is assured by an electrode 
structure, typically one ring electrode and two endcaps, that creates a static electric 
quadrupole field at the centre of the trap. Radial (xy) confinement is obtained by a static 
magnetic field B along the trap axis (z). Once trapped, a charged particle oscillates 
between the two endcaps with 'axial frequency' uj z . In the radial plane there are two 
motions with distinct frequencies: a fast motion with 'modified cyclotron frequency' u + , 
and a slow motion with 'magnetron frequency' The latter represents a slow drift 
about the centre of the trap. The sum of these two frequencies equals the 'cyclotron 
frequency' uj c = qB/m, where q is the charge of the single particle and m its mass. 

The electric field inside the Penning trap is described by the quadrupole potential 
U, which is given by 



where U is the trapping potential. The constant Rq = r$ + 2zq takes the geometry of 
the electrodes into account, where zq and are the distances from the trap centre to 
the endcap and the ring, respectively. A convenient choice for the trap dimensions is 
Tq = 2Zq, so that the trap depth is Uq/2 [7j. 




(1) 
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In the axial or ^-direction, the equation of motion is 



mz = —q 



dU 

dz 



Rl 



From equation (2) the axial frequency can be directly obtained, i. e. 
2 MU 



(2) 



(3) 



In the radial or xy-plane, the equations of motion are coupled by the magnetic field, 



i.e. 



mx = q 



my = q 



^-x + ijB 



R 



y — xB 



(4) 



It can be shown that the modified cyclotron frequency u + and the magnetron 
frequency oj- are given by |H] 



(5) 



A more detailed discussion of single particle motion in an ideal Penning trap is given in 

13 El El- 



3. Image charge effects 

The above description only holds true for the idealised case where there are no 
interactions between the particle and its environment. Therefore we need to consider 
the interaction between a single charged particle and the 'image charges' induced in 
the trap electrodes (this section and ^Oj). When there are more particles in the trap, 
also particle-particle interactions or 'space charge effects' have to be considered (section 
4). In order to understand the formation of image charges we will study two different 
geometries: a charge between two parallel conducting plates, and a charge between 
two hollow conducting spheres. It is worth pointing out that electronic detection and 
resistive cooling are based on the detection of image charges via a resonant RLC-circuit 
attached to the trap electrodes (sections 5 and 6). 



3.1. Charge between two parallel conducting plates 

The interaction between a charged particle and a single conducting plate can be 
described via the 'method of images' [H1II2]- A particle with real charge q, positioned 
at a distance z in front of the surface of an infinite conducting plate, induces a surface 
charge density a(r, z) which can be calculated by Gauss' law and is given by 

q z 
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Figure 1. A real charge between two parallel conducting plates and its induced 
images. The arrows illustrate the successively induced images qi at distances from 
the plates, R is the plate radius and d the plate separation. 



+ y 2 and (x, y) is the position on the surface. On an infinite plate the 



where r 2 = x' 
charge induced within a radius R is given by [13 J 



R 



a(r, z)2nrdr = —q 1 



(7) 



VR 2 + z 2 , 

such that for R — > oo the induced charge is g« = —q. 

A charged particle located between two parallel conducting plates induces a series 
of images. This is schematically shown in figure ^ by the curved arrows. Upon a 
displacement z, the real charge q induces an image charge q\ in the left plate at z\. 
This, in turn, induces an image q R in the right plate at z R , and so on. A second 
sequence of images also appears, starting with gf induced in the right plate. For a 
positively charged particle (q > 0), positive (even) images are located at even multiples 
of the plate separation d, i.e. the distances ±2d, ±4d, .... Negative (odd) images are 
located at odd multiples of d, i.e. the distances ±d, ±3d, .... The top dashed lines 
indicate the positions of the image charges for z = 0, the lower ones (and the small 
black arrows) for z > 0. 

From equation (7) it is clear that the image charge depends on the radius R of 
the plates, which implies that for finite R the image q^ is 'imperfect'. For an image 
% = 1,2,... induced in the left plate, the position z\ and the corresponding image 
charge qf are given by 

d 



(2« - 1; 



and 



Qi 



1 - 



R 2 + (z^) 2 / 



(8) 



where q$ = q$ = q. The expressions for the images induced in the right plate are 
similar, but with the sign of z reversed. Here we have made the assumption that the 
charge induced in a finite plate of radius R is the same as that induced within a radius 
R of an infinite plate. When R is small compared to d, the images are imperfect and 
the series terminates rapidly. The induced charge difference Aq between the two plates 
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Figure 2. Induced charge difference Aq plotted versus the position of the real charge 
q centered between: a) two parallel conducting plates (dashed line) , and b) two hollow 
conducting spheres (solid line). The black dots are calculated using equation (9) for 



z > 0. 



is defined as Aq = q L — q R , where q L and q R are the sums of the charges induced in 
the left and in the right plate, respectively. When the particle is located close to the 
right plate, Aq is positive, and vice versa. For a finite plate radius R < d, it suffices to 
consider only 3 images induced in both plates, as shown in figure HJ to obtain a good 
description of Aq. For such a simple case, it can be shown that Aq is given by a power 
series with only odd powers of z, i.e. 



The expressions for 7^ show that mainly the odd images contribute to Aq and that they 
lead to (higher order) odd powers of z. For increasing R, the variation becomes close to 
linear, which can be verified by extending the series to infinity. In that case the higher 
order odd terms vanish almost completely, leaving Aq « jiz. 

In figure |2l the calculated dependence of Aq is plotted as a function of the position 
z of the real charge q = 1 between two parallel conducting plates (dashed line) with 
radius R = 0.5 mm, which are separated by a distance d = 2.0 mm. Figure El clearly 
shows that for small z, Aq is almost linear with z, whereas for larger z, Aq varies with 
odd powers of z. For fixed R, when d w R, Aq is very non-linear at large z, due to 
strong higher order terms. When d > R, Aq is rather linear, but the induced charge 
difference is small. For infinitely large plates R — > 00, the response is very linear. Below 
we will show that higher order odd terms not only occur for planar electrodes, but also 
for curved electrodes (solid line in figure |2J). 

3.2. Charge between two hollow conducting spheres 

In a conventional Penning trap geometry, the endcap electrodes are better approximated 
by two hollow conducting spheres, ignoring the presence of the ring electrode. The case 



Aq « 7! z + 73Z 3 + 7 5 ^ 5 + . . . 



(9) 




Figure 3. a) Schematic of a hollow conducting sphere used to obtain the relationship 
between the real charge q and its image Q. b) A real charge q between two hollow 
conducting spheres and the corresponding first image charges q^, qi and q^ . 



of a charge outside a single hollow conducting sphere can easily be solved (see e.g. [T3]). 
As illustrated in figure Et), the real charge q, located outside the sphere at X, induces 
an image charge Q located at r inside the sphere. 

The potential at a point P on the surface of the sphere, is given by ^Hj 



v 1 



1 + i 

a b. 



(10) 



47re 

where the distances 

a 2 = R 2 + X 2 + 2XRcos(a) 

b 2 = R 2 + r 2 + 2rR cos(a) (11) 

are given by the cosine rule, and e is the permittivity of free space. After satisfying 
the boundary condition that the potential at any point P on the grounded conducting 
surface is zero, the real charge and its image are related as follows 

r = R 2 /X 

q= -QR/X (12) 

Now two spheres A and B with radius R and separated by 2D, as depicted in figure 
Eb); are considered. Upon a displacement z, the real charge q induces an image charge 
qf in sphere B at zf . This, in turn, induces an image q^ at z^ in sphere A, and so on. 
A second sequence of images also appears, starting with q± induced in sphere A. For 
the image i (integer) induced in sphere A, the position zf and the corresponding image 
charge qf are given by 

A 7-, 1 

Za = D — — j— r- and 



(D 
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-Qt-iR 



(13) 



where Zq = z A = z and q B = q A = q. The expressions for the images induced in 
sphere B are similar, but with the sign of z reversed. In analogy to the case of parallel 
conducting plates, the induced charge difference is defined as Aq = q A — q B , where q A 
and q B are the sums of the charges induced in sphere A and sphere B, respectively. Also 
in this case, for a finite sphere radius R, it suffices to consider only 3 images in both 
spheres to obtain a good approximation to Aq. Again, for such a simple case, it can be 
shown that Aq is given by a power series with only odd powers of z. 

In figure 121 Aq is plotted as a function of the position z of the real charge q — 1 
between two hollow conducting spheres (solid line) with radius R = 0.5 mm, which are 
separated by a distance 2D = 3.0 mm, giving a gap of d = 2.0 mm. The black dots 
are calculated using equation (9) for the 2 spheres case, taking only 3 images in both 
spheres into account. From equation (9) in figure |2 it can clearly be seen that the higher 
order odd terms are non-negligible for larger values of z. This analysis shows that also 
here, as in the case of 2 plates, mainly the odd images contribute to Aq and that they 
again lead to higher order odd powers of z. When D ph R, Aq is very non-linear at large 
z, due to strong higher order terms, but when D > R, Aq is rather linear. 

4. Frequency shifts due to image charge and space charge effects 

The image charge and space charge effects change the effective trapping potential and 
therefore lead to frequency changes. In the following, starting from the potentials at 
the trap centre, we derive analytic expressions for these frequency shifts and compare 
their predictions with experimental results obtained in measurements of clouds of ions 
confined in a Penning trap. In the case of the space charge effects we have also performed 
simulations, which confirm and support our analysis and are in line with experimental 
results. 

4-1- Image charge shift 

As can be seen in figure ^ the effective distance between the particle and the positive 
(even) images does not change as it moves between the plates. Therefore, only negative 
(odd) images exert a net force on the particle. For a single confined charged particle, 
the total force is thus given by the sum of the force due to the quadrupole potential and 
the image force. From figure Q it can be seen that the image force is given by 



where the summation is over odd images i only. Wineland et al. ^3] stated that 
the image force is given by _F* m « — 8q 2 z/(4neod 3 ), under the assumption that the 
displacement z is much smaller than the separation between the plates d. Our analysis 




(14) 
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shows that equation (14) is better approximated, even for small z, by including the 
higher order terms, i.e. 

Fr = Z^B(L + I + I + (15) 
47re Lrf 3 VI 3 3 3 5 3 )\ v ' 

Obviously, for larger amplitudes, i.e. z ~ d/2, this approximation is no longer accurate. 

In general, equation (14) is much better approximated, and over a larger range of z- 

values, when represented by a power series expansion. Ignoring the small higher order 

terms in equation (15), the image force is given by 

(16) 



-q 2 



8z 64z 3 
d 3 + ~d^ + 



4vre 

In analogy with section 3, it can be shown that mainly odd images contribute to the 
image force and that they lead to (higher order) odd powers of the particle displacement 
z. However, since the higher order effects are small for harmonic motion (small z), we 
will here just consider the first order term in equation (16). The corresponding equation 
of motion in the axial direction is given by 

4qUa z | 2q 2 
Rq ire d 3 
The image charged shifted axial frequency is thus 



mz = —^-z + M z (17) 



CO 



112 



4qU 2q 2 



To first order, the shift itself is given by 

1 a 2 

Au z = J' z -u z ^ ^— (19) 

7reo mw z d 6 

where the higher order effects have been neglected. If there are more particles between 
the plates, the image force and thus the axial frequency shift scales with the number of 
particles N. A result similar to equation (19) is obtained when the particle is considered 
to move inside a hollow conducting sphere [SI • The image force affects the centre- 
of-mass(c.m.) motion of a particle cloud, because it is an induced external force acting 
on the cloud clS 8b whole. 

Due to its radial motions, i.e. the magnetron and modified cyclotron motions, the 
charged particle also induces image charges in the ring electrode, which affect the radial 
frequencies. In order to calculate this radial image force, an approximate treatment 
similar to that given in section 3.2 can be applied, where the ring electrode is a small 
central section of the surface of a sphere. 

As can be seen from figure Ek), the force between the charges q and Q is given by 

™m 1 QQ 



r ~47re (X-r) 2 
By substituting equations (12) into (20) one obtains 
_,_ -a 2 R 



(20) 
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which is the radial image force on the particle due to one sphere. However, now the 
real charge q is located inside the sphere and its image lies outside. As the ring is 
considered as a small central section of the sphere's surface, for this case there will also 
be an image force like equation (21), albeit lower in magnitude due to the fact that 
not all the field lines originating from the charge q are directed at the section. Or, in 
other words, the image charge that is formed is imperfect. A charged particle with, for 
example, a magnetron radius r will thus experience a shifted magnetron frequency due 
to the induced charge in the ring. 

For a positively charged particle, the radial image force is directed radially outward, 
following the direction of the electric quadrupole force. The corresponding shifted radial 
frequencies are given by 

^ = Y ± nY) -r*-^ (22) 

where 7 is related to the ratio of the surface areas of the ring and the sphere. 

The image charge shift is only significant for very small traps and in most cases 
represents much less than 1 % of the axial frequency. Measurements of the image 
charge shift dependence on the number of particles have been performed for 1 H + , 
2 H + , 3 He + , 3 He 2+ , and 12 C 4+ ions [3] and frequency shifts were observed even for 
small numbers of particles. The magnetron frequency was shifted upwards and the 
modified cyclotron frequency was shifted downwards, as expected from equation (22). 
The relative magnitudes of the observed shifts, about 20 mHz per ion per charge, are 
also in fair agreement with those predicted by equation (22). For the cyclotron frequency 
oj c = u> + + u;_ no shift has been observed, as expected. 



4-2. Space charge shift 

Similarly to the image charge effect, there is a potential due to the presence of space 
charge which adds to the trapping potential and thus changes the trapping frequencies. 
The space charge potential U' is the potential seen by a particle due to the presence of 
all other particles and can be derived in a manner similar to the quadrupole potential. 
As a simple example, it is assumed that the particle cloud can be described by a perfect 
sphere with radius R' and homogeneous charge density p = Nq/V, where N is the 
number of particles and V = 4ttR' 3 /3 the volume of the sphere. Laplace's law in 
spherical coordinates (rj, 9, (f>) applied to this particle cloud states that 

^+2^: = _p (23) 

drf rj Orj €q 

where spherical symmetry with respect to 9 and holds. 

Outside the sphere the potential can be described by that of a point charge, i. e. 

Nq pR' 3 , . 

U'(r] > R') = —2- = % — 24 
Ane ri 3e n 
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Inside the sphere the potential should vary with if in order to satisfy Laplace's equation. 
The final form of the potential is 

U' = ~^(3R' 2 - rr 2 ) (25) 
6e 

For a single particle, the equation of motion in the ^-direction, with the substitution 
T) —> z, becomes 

mz = - q -(U + U) = -^ f z + -z (26) 

Both terms on the right-hand side of equation (26) vary linearly with z, i.e. equation 
(26) has the same form as equations (2) and (17). Therefore the space charge shifted 
axial frequency is given by (see e.g. [IE]) 

^ = ^-^- = ^ z -4 (27) 



where u 2 = pq/(me ) is the so-called 'plasma frequency' (see e.g. [17J. The factor 3, 
which appears in the denominator of the space charge term in equation (27), accounts 
for the geometry of the particle cloud. In a more general case, when the particle cloud 
is not a sphere but rather an ellipsoid, this constant changes. If the cloud shape is 
'prolate' the factor is smaller than 3, if the shape is 'oblate' it is larger PE]. To first 
order, the absolute frequency shift is given by 

2 

Ato z =u' z -u z & - Uq (28) 
bmeoUz 

where we have used the particle number density n = N/V = p/q. The relative frequency 
shift 

AlJz_ _ u)' z -u z _ nqR 2 
lo z u z 2Ae U 

depends linearly on the charge density p = nq and thus, for a given cloud size, on the 
total number of ions N = pV. 

The space charge limit is reached for \U'\ > \U\, i.e. when the space charge potential 
compensates the trapping potential. In this case there is no potential minimum and 
confinement is lost. 

If the space charge density is not homogeneous, but has a Maxwell-Boltzmann 
distribution, the space charge shift depends on the average distance of the particle 
from the cloud centre [TH]. This leads to a distribution of shifted axial frequencies. 
However, the c.m.-motion of a particle cloud is not affected by intra-cloud interactions 
and thus remains unshifted by space charge effects 14J . A similar broadening of the axial 
frequency distribution occurs when the trapping potential is anharmonic. In that case, 
the axial frequency depends on the amplitude of the motion such that a distribution of 
particle energies automatically leads to a distribution of axial frequencies. To correct for 
such anharmonicities, most Penning traps have compensation electrodes (see e.g. |19j). 

The space charge potential U' also affects the radial motions and thus the radial 
frequencies of the particle. Due to the difference in signs between the z- and r- 
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dependences in equation (1), the space charge term now has the opposite sign. This can 
be shown via the substitution U — > U' + U in equation (2). The space charge shifted 
radial frequencies are given by 



Thus, there is a downwards shift of the modified cyclotron frequency u + and an upwards 
shift of the magnetron frequency u;_ with increasing charge density p. 

Measurements of the magnetron frequency of an ion cloud in a Penning trap with 
tq = 5 mm, zq = 3.5 mm, Uq ~ 6 V and B = 1 T were performed by Dholakia et 
al. jl] through laser cooling and laser spectroscopic techniques. For a single ion they 
yield a magnetron frequency of ~ 40 kHz, while for clouds with ion number densities of 
the order of 10 7 cm -3 , the corresponding value is shifted to ~ 60 kHz. This is in fair 
accordance with the behaviour expected from equation (30). 

4-3. Simulations of the space charge shift 

In an effort to check the validity of the results of the calculation, the influence of space 
charge has been studied using SIMION, which allows for charged particle tracking in 
static electric and magnetic fields created in a user-defined geometry. The simulations 
have been performed for hydrogen-like lead ions EH] = 207 u and q = 81) confined 
in a cylindrical open endcap Penning trap f!9| . 

The ion cloud is simulated by a cubic array, typically of size 5x5x5, in which a 
charge of magnitude e' is assigned to each array point. So, for a cloud of N = 125 
charged particles e' = q x e (e = 1.6 x 10~ 19 C), whereas for a cloud of N = 125000 
particles, e' = q x 1000 x e. The particle-particle interactions are simulated by SIMION 
using the grouped Coulomb repulsion e'. Within this approximation, small clouds are 
clearly better described than large ones. It may be speculated that a simulation where 
N particles are simulated by N array points would result in the frequency shift expected 
from equation (29). This, however, is not within our computational capabilities since 
the computation time scales with iV 2 . Once a simulation produced several periods of 
oscillation (for all particles), the simulation was terminated and the particle frequencies 
were extracted from the data by means of Fourier analysis. 

The results of our simulations are shown in figure where the relative axial 
frequency shift is plotted versus the number of particles in the cloud. The parameters 
were: z = 15 mm, Uq = 1000 V, B = 1 T, JV = 125000 (5x5x5 array) and n = 4 x 10 6 
m~ 3 . The simulations show a linear relationship between the relative axial frequency 
shift and the number of particles. Also shown is a plot of equation (29), which behaves 
linearly for small n (oc N) as expected, but clearly deviates from the results of our 
simulations for larger numbers of particles. This discrepancy is attributed to the limited 
number of array points used in our simulations, the cubic lattice geometry (rather than 
spherical), and the method used to simulate high ion densities (i.e. via e'). However, 




(30) 
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Figure 4. Simulated relative axial frequency shift —Aui z /ui z plotted versus the number 
of particles N. Equation (29) is also shown for comparison. The inset shows the 
dependence of the shift on the distance from the cloud centre. 



for large N also the approximation used for the derivation of equation (29) is no longer 
accurate. 

The c.m.-motion of the cloud is found to be invariant under the influence of 
interaction between the particles, as expected. However, as shown by the inset in 
figure EJ the relative frequency shift is larger for particles that are initially positioned 
close the cloud centre than for particles that are further out. This is attributed to the 
fact that the amount of space charge surrounding the particle drops with increasing 
distance from the cloud centre. 



5. Electronic detection of the axial motion 



In order to experimentally determine the axial frequency, a resonant RLC-circuit with 
a high quality factor (see section 6) is connected to the endcaps. A Fourier analysis of 
the induced current / in the RLC-circuit, picked up by a sensitive amplifier, then yields 
uj z . For the case of a charge between two parallel conducting plates, the current flow at 
any time is 

dAq dAqdz dAq 

1 = ~dT = ~dVdi = Vz ^z~ {d} 

where Aq is the induced charge difference between the two endcaps. For parallel plates 
with large radii, Aq varies linearly with z, i.e. Aq = jiz. However, in many other 
cases, as we have shown above in equation (9) and figure El there can be higher order 
odd terms as well. For example, if Aq = 71 z + 73Z 3 , it can be shown that the electronic 
circuit also picks up the higher order term. For simple harmonic motion, of the form 
z = zo cos(uj z t), one obtains 

Aq = 7xz cos(u z t) (32) 
+ ^-2- [3 cos (oo z t) + cos (3oo z t)] 
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The cubic term therefore induces signals at u z and 3uo z . 

Consider a single particle with mass m and charge q which is moving inside a 
particle cloud of mass M and charge Q. In this case only the c.m.-motion of the cloud 
will normally be detected. If the particles are all identical, then q/m = Q/M. If the 
position of q is z and the position of Q is Z, then Z = —mz/M = —qz/Q in order for 
the cm. to be stationary. This would be the case for small oscillations of the particle 
about an equilibrium position, i.e. the normal modes of the system. The induced charge 
difference, with q and Q moving such that the cm. is stationary, is 

Ag = i 3 z 3 - ^ « l3 z 3 (33) 

Again the signal is induced by the cubic term, which gives rise to signals at uj z and 3u z . 
If the oscillation frequency of q in the environment Q is reduced from the single-particle 
value, then this shifted frequency will be seen in the signal around u> z , so long as the 
term in 73 is present. 

If the amplitudes of motion are larger than the interparticle spacing, the normal 
mode description no longer applies. In this case the motion of an individual particle does 
not necessarily give rise to a compensating motion of the rest of the cloud. Therefore, 
the individual oscillation frequencies of all particles are observed directly in the detected 
signal i-e. via the non-linear term in equation (32). 



E(t) = £ exp(-t/r) and r = — ^f- (34) 



6. Resistive cooling of the axial motion 

For a single ion, the current I running through the resonant RLC-circuit results in an 
exponential decrease of the particle's kinetic energy, which is given by ^1] 

m(2z ) 2 
Rq 2 

where R is the resistance of the RLC-circuit. In resonance, the impedance of an RLC- 
circuit is real and acts as an ohmic resistor with resistance R = QLu z , where Q is 
the quality factor of the circuit and L the inductance. The bandwidth is given by 
Au z = v z jQ and for typical values [y z ~ 1MHz, Q ~ 1000) the bandwidth is ~ 1 kHz. 
The resistive cooling time constant, expressed in terms of the circuit components, is 
given by 

4(z 0y /m) 3 



LQVrVo (35) 

The cooling process finally leads to an equilibrium between the particle's kinetic energy 
and the temperature of the heat bath. However, the noise temperature in the electronics 
may exceed the physical temperature of the RLC-circuit, leading to somewhat higher 
temperatures of the particles Deviations from the exponential cooling behaviour 
may occur when the dissipated power is not proportional to the kinetic energy of the 
particles. For example, if the trap is anharmonic, uj z depends on the particle's kinetic 
energy and may move out of resonance with the RLC-circuit 0111]. 
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For coherent c.m.-motion of a particle cloud with N particles, the total induced 
current is given by I tot = NI, and correspondingly this motion is strongly cooled. Apart 
from the coherent c.m.-motion, there can be non-coherent motions of ions within the 
cloud. As described above, these individual ion motions give rise to a resultant induced 
signal in the RLC-circuit \17\ 122] . Since the signals from the non-coherent motions 
cancel to first order, the corresponding cooling is weaker and therefore the cooling times 
are longer |24j . 

As discussed in section 4.2, space charge effects generally lead to a shifted and 
broadened distribution of axial frequencies of trapped particles. Resistive cooling is 
only effective within the bandwidth of the external circuit, which may be less than the 
range of frequencies expected from the space charge effects. In large ion clouds the axial 
frequency distribution may be considerably broadened by space charge effects (see e.g. 
inset figure EJ). In that case, it is expected that the inter-particle collisions continuously 
thermalise the cloud by redistributing energy over the cloud's axial frequencies so 
that the resistive cooling process continues, but at a reduced rate limited by the 
thermalisation process. The inter-particle collision rate for a single component, trapped 
particle cloud is characterised by k oc g 4 m _1 / 2 nT -3 / 2 where q is the particle's charge, 
m its mass, n the particle number density and T the intrinsic temperature of the ion 
cloud Therefore, thermalisation is most efficient for light, highly charged particles 
in a dense, cold ion cloud. It may be expected that, for an efficient thermalisation 
process, the cooling time of a particle cloud scales with the bandwidth of the RLC- 
circuit. Therefore, in practical applications, a compromise has to be found between 
a sufficiently small cooling time constant r oc 1/Q and a sufficiently large bandwidth 
8v z oc 1/Q. 

7. Conclusions 

We have presented a thorough analysis of the formation of image charges. Such a study 
is relevant because the images formed in the trap electrodes of a Penning trap will be 
picked up by an electronic detection circuit attached to the trap electrodes. Our analysis 
was done for the case of two parallel conducting plates and for the case of two hollow 
conducting spheres. It was found that in both systems mainly the negative (odd) images 
contribute to the induced charge difference Aq between the two electrodes. These images 
lead to (higher order) odd dependences of Aq on the particle's position z between the 
electrodes and its oscillation frequency uj z . We therefore concluded that an electronic 
detection system can, in principle, pick up higher order frequencies of the particle motion 
as well. Such frequencies can, for example, be induced by individual particle motions at 
a frequency 3u z , rather than exclusively at the common centre-of-mass motion at the 
frequency uj z . 

We have also derived analytic expressions for the shifted frequencies of trapped 
particles, which are due to image charge effects and space charge effects. Both effects 
lead to a downwards shift of the axial frequency uj z and the modified cyclotron frequency 
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lo + , but to an upwards shift of the magnetron frequency These shifts depend strongly 
on the charge density of the confined particle cloud and on geometrical trap parameters. 
Our simulations qualitatively confirmed the downwards shift of the axial frequency. In 
addition, they also showed that in cold ion clouds, with motional amplitudes smaller 
than the cloud dimensions, the frequency shift depends on the particle's average distance 
from the cloud centre. 

Resistive cooling of a trapped particle cloud depends also on the image charges 
induced in trap electrodes, and on the energy dissipation in an external heat bath 
formed by a frequency resonant RLC-circuit. Due to the limited bandwidth of such 
a circuit, space charge induced shifts and broadening of the axial trapping frequency 
distribution need to be taken into account. Furthermore, to first order only the centre- 
of-mass mode induces a signal in the circuit. Since the cooling is effective only within 
the bandwidth of the circuit, resistive cooling of large clouds, where considerable space 
charge effects occur, may differ significantly from the single particle case and may lead 
to much longer cooling times. 
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